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The overall goal of this project, was to nnderst ami in more detail how a pattern-burning system can 
adjust its spacing. ‘'PaUern-lorming systems/' in this context, are noneqnilihr iiun contina whose slate is 
determined hy some experimentally adjust a hie control parameter. Below some critical value of the control 
parameter, the system remains in a simple steady state. Above critical, this state becomes unstable, and the 
system then has available to it a range of linearly stable, spatially periodic st eady stat es, each characterized 
by a spacing which can lie anywhere within some band of values, t hese systems include inteiTacial systems 
like directional solidification, where the solidificat ion front is planar when the ratio of growth velocity to 
thermal gradient is below its critical value, but takes on a cellular shape above critical. They also include 
systems without interfaces, such as Beuard convection, where it is the fluid • locity held which changes from 
zero to something spatially periodic, as the control parameter is increased through its critical \alue. 

The basic question, to be addressed was that of how the system chooses one of its myriad possible 
spaeings win u t he control parameter is above critical, and in part icular the role of noise in this selection 
process. Previous work on explosive crystallization had suggested that one spacing in the range should 
be preferred, in the sense that, weak noise should eventually drive the system to that, spacing. 1 hat. work 
had also suggested a heuristic argument for identifying (lie preferred spacing, the project had three main 
objectives: to understand in more detail how a pattern-forming syst cm can adjust its spacing; to invest igate 
how noise drives a system to its preferred spacing*, and to extend the heuristic argument for a preferred 
spacing in explosive cryst allization to other pattern-forming systems. 

The project began with the first of these objectives. With an undergraduate student, .led Overmami. 
I examined the process of pattern adjustment in the Swift- llohenberg equation, which is a standard model 
of interfacial pattern formation, and some of its generalizations. 1 hose equat ions have linearly stable steady 
states of the form Tsin(A\r) for a range of wave number k. We developed numerical algorithms to find the 
“saddle” states, the ones which are intermediate between the steady stale having N waves on a certain 
interval and the one having N T 1 waves. After some exploration, we came to realize that these saddle states 
are essentially phase perturbations of the steady st ates, in which two ol t he N waves are pulled fart inn* apart 
than the average spacing, with several neighboring waves pushed very slightly closer together. l)epriidmg 
on the detailed shape of the interface in the gap between t he separated waves, a new wave may or may not 
form there, resulting in a return to the Af-vvave state or an adjustment to the A -|- 1-wave state* respectively. 
The relevance of this work to the overall project will be discussed below. 

To address the second objective, I investigated the fokker- Planck equation, which governs the probability 
distribution for states of a system subjected to noise. Let tin* front shape he given by its Fourier spectrum, 
with the component of wave number /■ having amplitude A^. I he evolution of the front shape is specified 
by equations for the t ime evolution of the 

A k =l\[Al 0) 


The notation A is a shorthand for the collection of all the /!/., and the dot represents a time derivative. The 
possible steady state front shapes are those for which all /q. vanish, for each Ay we now add an indepemh ut 
Gaussian noise source with mean zero and variance ?/“. I his is the appropriate type of noise to add in order 
to test which steady slate of ( l) is preferred by the dynamics, since by giving each wave number equal power, 
the noise itself does not favor any wave number. The front shape A then becomes a random variable', and we 
seek its probability distribution T [ A ] . AfI.er an initial transient, I his is given hy the stationary fokkei I laiuk 


equation, 
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with summation over repeated indices implied. Note that (lie noise strength 
derivative in the equation, so that the weak-noise limit, ?/ — * 0, is singular, lo 


write the P in a WKB form, 

V[A]= /v [A] ex ]>( — IF [.■]]/?/) ( 1 T 0{i})) 


i] multi ply ios the high 
solve (2) in 1 his limit , 
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Substituting this into (2) and expanding in powers of ?/ yields equations for W and /v . From {V, we sec that 
the most probable state of the system in the weak- noise limit is that for w inch il is a minimum. I bus the 
correct way to determine the preferred stale of the system is to compute 11 lor each steady state, time find 
the one for which it has the lowest value. 

An important special case arises when the deterministic equations (1) have a orathnit shuet itu , i.e.. 
when we have /').[.!] — — iW/DAt for some scalar function In lliis case, we have I he exact wlulinii 
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p — f*xp(— 2<l>[/1 j / r/ ) . Although most of i\\ n systems wr ;u<' into rostral in do not. h;»vn this stnif(un\ systems 
whirl] do provide useful rlirrks on the cab/ubil ions. It is easy to see that, (ho (U Icrwinisfir rvulufinn (I) 
makes $ decrease with time until the system reaches a steady state. 'I Inis the linearly stal.de steady stales 
are local minima of d>. We found that IT plays almost the same role for systems lacking a gradient d.mdiiw: 
the deterministic evolution never increases IF. In part icular, linearly stable steady st-ates of t he tleienninist-ic 
dynamics are local minima of IK. 

These general considerations show that we need to calculate H [A\ appearing in (•>). 1 he equation b>i 


'Ml' /l/MK 
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This is the Hamilton-.Jncohi equation at zero energy for the Hamiltonian V\f)\ T] \f>\ ? ‘/2 I /v,U. or the 
Lagrangian £[<IA/(H , A] — (L/2)(/i*. — /q-)(/U ~ f'fcb u 'ilh II Indiig the action. 



Fo calculate the difference in IK between two states, we should find t he least. -act, ion path conneci-iiig those 
points and integrate along it. according to (•)). Least-action paths can be found by integrating the Laginnge 
equations in configuration space or the Hamilton equations in phase space. However, since C is manifestly 
positive, W must increase along the path. Hut a stable steady state is a local minimum of IK, so least-action 
paths must all come out of such a state, dims there cannot be a path from one stable steady state to another. 
To compare two steady states, we must integrate paths running from each steady state to some common 
point. 

To clarify what is involved in these calculations, we have written numerical code for a simple case with 
only two degrees of freedom and two attractors, t he deb monistic equations are 

x — x — x 3 — bxu J \ 

" „ (b) 

if -z rtf / — cy* — dxij. 


These equations have a gradient, struct tire only for the special case b — <!. I hey arc' obtained by taking a 
typical pattern-forming system in the weakly nonlinear regime (control parameter slightly above critical) and 
keeping only two Fourier modes, whose amplitudes are denoted here by .r and if. All details ol the particular 
model and the choice of which wave numbers to keep are embedded in the constants h, c, d, ami e. Notr 
that these equations arc not meant to he a serious model of any real pattern-forming system. Hath* i. we 
have studied them to clarify the issues that will arise when doing calculations for real systems, calculations 
which will need to be carried out in high-dimensional, rather than two-dimensional, configuration spaces. 
The model in (6) has two attractors, one on the ;?• axis and or- on Hie y axis, whose basins are separated 
by a separatrix which contains a repel lor at the origin and a saddle point.. Wo found that, theie an' paths 
satisfying the Lagrange equations which run from each attractor to the origm. but (or b *f one of those 
always crosses a caustic, and so is not a minimum of the action, thus in older to compute the 11 \ulms 
at the two attractors, we must find the least-action paths running from the attractors to tim saddle, ami 
so compare W at each attractor with W at the saddle. Note that a by-product of t hese calculations is the 
relative values of IK between Hie attractors and the saddle. This turns out to give the exponential noise 
dependence of the rate of escape from each attractor, and so an estimate of the lifetimes n| the linearly stable 
steady states. 

The next step, which is currently under way, is to lake what wo have learned from our anahtio.il 
calculations, our low-dimensional numerical calculations, and our search for saddle st all's, and devise pi a. I d al 
numerical methods for comparing the IV’ values of attractors in much higher dimensional systems. I he 
attracting states are. easy to calculate, since the time evolution of nearby states converge* to them. Wo 
have developed methods of calculating the saddle stales. VVliat remains is to develop ellinent methods (or 
finding the least -action path from an attractor to its neighboring saddles. Two types of method are possible. 
Relaxation methods start with a trial path from the attractor to the saddle, then vary the path m order 
to minimize the action. With n points along the path in a (/ dimensional space, this is a minimization over 
vd variables. The alternative is a shooting method, starting at the attractor, choosing an initial direction 
for the path, and integinting the Hamilton or Lagrange equations, then varying the initial direction m an 
attempt to hit the saddle. This is a. minimization over only <1 - l variables, but each trial takes much mote 
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computation than moving a point along a trial path in the relaxation m'hhnii. 1 !i p relaxation niHliml uses 
the fact that the solution path must go to the .saddle, hut not the fart that if. iimsl also satisfy the hngrange 
equations, while tin' shooting method uses the latter fad- hut not. tin - former, this is one of two dived ions 
this research has lv-nn f aking since the conclusion of the contract. 

Clearly tie calculations which give the correct answer for the preferred spacing are diflirull and involved. 
It would be much better if there were some vastly simpler way of finding the preferred spacing directly, 
without needing to calculate a variety of intermediate results. I his is where the third objective u'liio m. 
A candidate for a simpler n; 1 aion arose in my earlier work on explosive crystallization; briefly, it is Mu' 
conjecture that the preferred spacing is the one for which the growth rate of the foiirier amplitude, rM*/f 7. 
regarded as a function of Ay has a maximum at the same wave number A lor which Ah i I sol f has a maximum. 
'Pile third objective was to explore the relation between this criterion an 1 the coned, one, in an attempt 
to determine whether and under what circumstances the two give the same result, i.o.. when the h< urisl ic 
argument is correct. 

We have been able (n show that, the heuristic criterion does indeed produce the correct preferred wave 
number when the deterministic equations have a gradient structure. The argument runs as follows. Suppose 
the free energy functional T[/l] is minimized among functions of the lorm ,4 ~ M^(h-tf) with a Jirr*l wave 
number (j by M — /\/q (//). Denote the corresponding value of T by 

>('/)*(■ ~ 7 )]- ( ‘ ) 


The preferred wave number will be the q value for which is a minimum, do find this, we calculate 

<P win (q 4- 6q)\ at the preferred wave mmiher, this will differ from <!’„„•,, (</) •>>’ an amount, which is of order 
( d/) 2 . We find 
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Tlio functional derivative and its /. - < I«*r i v;vf ive ;i|>pi';iriiig on tin' right side ol tin' final equation ate evaluated 
at. Ai- = M 0 {q)Hk ~ <i) and k -- </■ Thus the second (.crni on the right vanishes, since it. is just, the right si do 
of the deterministic equations, /q.[,l], evaluated at. a fixed point. In order tor q to bo the minimum of 
the third term on the right must then vanish also; in other words. i>I\-/Ok, viewed ns a Inaction of k. must 
have a. peak at k ~ q. 1 his is precisely what, the heuristic argument claims. 

It. is difficult, to see how this argument can Ire generalized to systems which do not have a. grad Hit- 
structure. Indeed, general analytical results of any kind are difficult to come by in the non- gradient case. I 
am currently trying to work out the argument when the deterministic dynamics is dose' (o a gradient system. 


that is when we have 
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where e is a small parameter. Adding the perturhation f»Y generally shills the positions of tie' fixed pu 
the dynamics, the deterministic trajectories, and the action-minimizing paths. Inking all this into ar 
we find that the difference in action between an attractor /f ,s) and a. saddle .•i'"' is given, to first ordr 

hy • r - 

1V’[/V ( '')] - If = 2'l , [,l"‘ ,n) ] - 2'f’[.'f ( ' ,n '] -2> r.Y[,l] .//l*, 

where and f * n! are. the mi; vtnrbcd positions of 1 lie saddle and attractor, and I lie line int> 

over the time-reversed, unperturbed deterministic trajectory which runs from the latter to the forme 
next, step, which is in progress, is to sock a relation between this result and the picdictiou e>l I hi' he 
criterion. 
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